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Abstract 

In this paper, we prove a Donsker type approximation theorem for the Rosenblatt 
process, which is a selfsimilar stochastic process exhibiting long range dependence. By 
using numerical results and simulated data, we show that this approximation performs 
very well. We use this result to construct a binary market model driven by this process 
and we show that the model admits arbitrage opportunities. 
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1 Introduction 

Long range dependence stochastic processes have been intensively used as models for different 
physical phenomena. First, these properties appeared in empirical studies in areas like hydrol- 
ogy and geophysics; more recently, they appeared to play an important role in network traffic 
analysis and telecommunications. As a consequence, efficient mathematical models based on 
long range dependence (or long memory) processes have been proposed in these directions. 

The notions of long range dependence and selfsimilarity have also been considered 
in mathematical finance. An excellent survey on the different aspects of the appearance of 
the long range dependence in practice is the paper [3]. The debate on the presence of long 



memory in stock prices is actually not new. The idea that asset returns could exhibit long 
range dependence comes from Mandelbrot ([llj) and then observed in several empirical studies. 
We refer, among others, to [26j for concrete examples and for interesting comments on this 
question. We also mention that some other authors rejected the idea of the presence of long 
memory in asset returns (see e.g. [TO]). 

A rather general opinion is that long range dependence in financial models is strongly 
related to the presence of arbitrage. For example, in the case of market models driven by the 
fractional Brownian motion, this has been explicitly shown by Rogers (|17j) or Sottinen (|22j). 
In special situations, for example under transaction costs, arbitrage could be eliminated (see [7] 
or [20]). Different approaches, based on Wick-Ito calculus, have been developed in e.g. [2], [9]. 
Models driven by long range dependence processes others than the fractional Brownian motion 
have been, from the stochastic calculus point of view, less considered. This is actually one 
of the motivations of our work: we propose a binary market model driven by a non-Gaussian 
selfsimilar process (called the Rosenblatt process) which exhibits long range dependence and 
we show that this type of model admits arbitrage opportunities. It is actually known that in 
general arbitrage- free models imply the fact that the price process is a semimartingale ([5]) 
but it has been recently proved in [4] that by modifying the class of admissible strategies, one 
can also consider arbitrage-free models driven by non-martingales. 

To construct our binary market model based on the Rosenblatt process we need a 
Donsker type theorem to approximate in law this process by some disturbed two-dimensional 
random walks; this results could be useful by itself. In fact, this theorem extends a result by 
Sottinen [22] and represents a variant of the so-called Non Central Limit Theorem proved in 
[6j and [2^. We mention that, since we are now in a non-Gaussian context, the proof of this 
result demands different techniques. 

Our paper is organized as follows. In Section 2 we describe the basic properties of 
the Rosenblatt process. Section 3 contains the proof of the Donsker theorem to approximate 
weakly, in the Skorohod topology, the Rosenblatt process by walks. In Section 4, we introduce 
our binary market model which is showed to converge to the Black and Scholes model with 
Rosenblatt noise. We show that the model admits arbitrage opportunities and we construct a 
such opportunity. 

Finally, in Appendix we present numerical results based on simulated data which show 
that the approximation method performs very well. Related numerical results can be found in 
[5 or [16]. 

2 Preliminaries 

The Rosenblatt process appears as a limit in the so called Non Central Limit Theorem (see [6] 
or [23])- We recall the general context. Consider (^n)nGZ a stationary Gaussian sequence with 
mean zero and covariance 1 such that its correlation function satisfies 




(1) 
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with k > 1 integer, H E 1) and L is a slowly varying function at infinity . Denote by 

Hm{x) the Hermite polynomial of degree m given by Hm{x) = {—l)"^e~ -^^e 2" . Let g he a 
function such that E(g(.^o)) = and E{g{^o)^) < 00. Suppose that g has Hermite rank equal 
to k; that is, if g admits the following expansion in Hermite polynomials 

g{x) = CjH,{x), Cj = |e (giCoHji^o))) , 

then 

k = min{j; Cj / 0}. 

Since E [5(^0)] = 0) we have k > 1. The Non Central Limit Theorem ([6], [24]) says that the 
sequence of stochastic processes 

[nt] 

converges as n — > 00, in the sense of finite dimensional distributions, to the process (called the 
Hermite process) 



zUt) = ciH, k) / {\\is- ' > dsdB{y,) . . . dB{yu), (2) 




where xj^ = max(x, 0) and the above integral is a multiple Wiener- Ito stochastic integral with 
respect to a Brownian motion B{y))y^M. (see p5]) . 

Let us list some basic properties of the Hermite processes. 

• it exhibits long-range dependence (the covariance function decay at a power function at 

zero -"Joseph effect") 

• it is i/-selfsimilar in the sense that for any c> 0, (Z^(rf)) =(^) {c"Z^{t)), where " =('^) " 

means equivalence of all finite dimensional distributions 

• it has stationary increments, that is, the joint distribution of {Z^{t + h) — Zfj{h),t G [0, T]) 

is independent of h > 0. 

• the covariance function is 

E(Z^(i)Z^(.)) = 1 {t^^ + .2^^ -\t- , s,tG [0, T] 

and consequently, for every s,t G [0, T] 



E 
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z^(t)-z^(.) =it-sr (3) 
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• the Hermite process is Holder continuous of order 5 < H 

• if /c > 2, then Zfj is non-Gaussian. 

When k = 1 the process given by ^ is nothing else that the fractional Brownian 
motion (fBm) with Hurst parameter H G (0, ^). For k > 2 the process is not Gaussian. If 
k = 2 then the process ^ is known as the Rosenblatt process (it has been actually named by 
M. Taqqu). 

We focus here our attention on the case k = 2. We will work with the representation 
of this processes as integral with respect to a Wiener process on a finite interval. Recall that 
the fBm with Hurst parameter H > ^ can be written as 

i?f = f K"{t,s)dW,, tG [0,r] (4) 

JO 

with {Wt-,t G [0,T]) a standard Wiener process and 

K"{t, s) = CHS-2-" j\u - s)"-^2u"-hdu (5) 
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where t > s and ch = (^ ^(2-2^^-1) / ' -^'^o™ ® obtain that for t > s, 

{t,s) = CH{-) (t-s)^-^. (6) 



Aa analogous representation for the Rosenblatt process has been given in [25j. We have 



zUt) ■.= Zt=^'^ d{H) f f 

Jo Jo 



dWy^dWy, (7) 



where {Wt,t G [0,T]) is a Brownian motion, 



H'^^^ (8) 

_ 1 

and d{H) = ( 2(2^-1) ) ^ • ^'^^^ -^^ > I implies H' > |. 

3 Convergence in law to the Rosenblatt process 

This part in consecrated to a Donsker invariance principle for the Rosenblatt process. From 
now on, we will consider the Rosenblatt process to be given by the formula ([7]). We will denote, 
for every t G [0, T] 
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and then ^ ^ 

Zt= [ [ F{t,yi,y2)dW{yi)dW{y2), t e [0,T]. (9) 
Jo Jo 

The kernel K^' (denoted simply by K in the sequel) is the standard kernel of the fractional 
Brownian motion ([5]). 

Let us first recall some known facts. Consider (^j)i>i a sequence of i.i.d random vari- 
ables with E(^j) = and IE(^^) = 1. The Donsker Invariance Principle says that the sequence 
of processes 

[nt] 

converges weakly, in the Skorohod topology, to a standard Brownian motion. Here [x] denotes 
the biggest integer smaller than x. 

This result has been extended in [22] to the fractional Brownian motion (see also [12j). 

Define 

K''(t,s) :=n r K{^^,u)du, n>l 

n 

and put 

B^= / K^{t,s)dW^ = y2n r KC^,s)ds^, n > 1. 
Jo ^ 'IT' Vn 

Then it has been proved in [22] that the disturbed random walk converges weakly to the 
fractional Brownian motion. 

From the above results and the representation ([9]) it is quite natural to define the 
following approximation for the Rosenblatt process 

Remark 1 We eliminate the diagonal "i=j" because the Rosenblatt process is defined as a 
double Wiener-Ito integral and as a consequence it has zero mean. When the diagonal i = j is 
included in the sum ill]) then the limit is in general a double Stratonovich integral (see or 



Proposition 1 The family of stochastic processes (^")tg[o,T] converges in the sense of finite 
dimensional distributions to the process (■^t)tg[o,T] ^W■ 

Proof: We will proof this result in several steps. 
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Step 1: Let us consider an arbitrary sequence of partitions of the interval [0, T] of the 

form 

TT : = Iq < ti < . . . < = 1 

with Ivr™'! ^ as m ^ oo. Define 

zr= E T^^\{j^jj(i^^^^)^-A^(^T)W (12) 

where we denoted by A™ = and by 

Ty(A™) = WtY^-WtYi^. 

Then it follows from |21j . Theorem 3.4, or [8] that for fixed t the sequence Z^"" converges in 
L^(r2) as |7r"*| — > to the multiple Wiener-Ito integral of F{t, •) with respect to the Brownian 
motion W t T 

[ [ Fit,u,v)dWudW, = Zt. 
Jo Jo 

Step 2: Secondly, define the process 

^r'"= i; ( f [ F{t,u,v)dvdn]w-{AT)W-{Af) (13) 

where is the random walk given by (jlOp . Then clearly, for fixed m, as n goes to oo, the 
finite dimensional distributions of Z^"^'^ converges to the finite dimensional distributions of 
Z'^™ (this comes from the weak convergence of to the Wiener process W). 



Step 3: We prove now that for every t S [0, T], the sequence Z^ converges in L^(r2) 
'j'" as m — > oo, where 

z>= ^ y y i *e[o,r]. (14) 



n n 



Consider the sequence 

F""'(i,n,^;) = E i^mLmi (/^/^-^(^'^'^)^^^^) lAr(^)lA™(^). 

Then F'^'"(t, •) converges to F(t, •) in L2([0,T]2) as m ^ oo (see [H], [S]). 
First note that Z^'" can be approximated in L?'{^1) as m — > oo by 

Z.' ' = > n / F (t,u,v)dvdu—=—= 
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Indeed, 

/ _|2 



[nt] / \ 2 



= E n^[l^l_^{F^"^it,u,v)-Fit,u,v))dvdu\ EiekM^f) 

k,l=l;k^l \ n n / 

[nt] k I 

< ^{F^^^ (t,u,v) - F{t,u,v)fdvdu 

k,l=l;kytl n n 

(F"" {t, u, v) - F{t, u, v)fdvdu 



and this clearly goes to zero as m ^ oo. 



It remains to observe that Z^'^ is equal to Zf"" for every t, m, n. We can write, if A 



denotes the Lebesque measure, 

[nt] k I 



zf = J2 ^ L 1 h 1 ^'"^^ 

> n^^^^^ / / Ia'^MIa'^ (v)dvdu 

k,t=l;k^l n n 

E "^|;JjA([t_i,i)nAr)A([i^.l)nA» 



[nt] 

X 

k,l=l;k^l 



^k 6 
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and on the other hand by using p2p and (jlO|) . one has 

1 



E 



|A"^||A™| 



m I Am 



and it is not difficult to see that Zj' 



and Z^''^ coincide. 



Step 4- At this point we conclude that the family of processes Z converges in the 
sense of finite dimensional distributions to the Rosenblatt process Zt- Let h he a function 
defined on W and consider si, . . . ,Sp G [0, T]. We will show that 



E(h{Z>,...,Z>))-E{hiZ,„...,Z,^)) 



converges to zero as n ^ oo. This can be bounded hy A + B + C where 



A 



E(/i(z,,,...,z,j)-E(Mzj™,...,z;; 



and 



B 



C 



E[h{z:^ '\...,zi 



E 



) • • • ) ^Sr, 



By Step 1 and Step 3 we have that for fixed c > the terms A and C are bounded (uniformly 
in n) by c when m is large enough, the term B tends to zero as n — > cxd from Step 2. 

Step 5: Clearly the family Z '" can be replaced by the family Z"(|lip because their 
difference goes to zero in Lp'{Q) as n ^ oo. ■ 

Next, we show the tightness. 
Proposition 2 The family given hy is tight. 

Proof: Let s < t, s,t £ [0, T]. It holds, since the kernel F{s,u,v) vanishes when u oi v are 
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bigger that s, 



< 



< 




, 'jit] \ [ [ns] 
n I \ n 



F[^^,u,v\-f{^^,u,v 



n 



n 



dvdu 



dvdu 



n \ n 



n I \ n 



dvdu 



T f-T 



JO 



F[^^,u,v]-f(^^,u,v 



n 



n 



dvdu 



[nt] [ns] 



n n 

Now the conclusion follows by using exactly the same arguments as in [22], end of the proof 
of Theorem 1. ■ 

The main result of this section is a consequence of Proposition [1] and Proposition [2l 

Theorem 1 The family of stochastic processes Ul\) converges weakly, in the Skorohod topol- 
ogy, to the Rosenblatt process. 

4 Binary market model based on the Rosenblatt process 

The binary market constitutes a financial model where the asset are traded at discrete times. 
In classical cases (for example when the driven process is the Wiener process) the binary model 
approximates the Black and Scholes model. 

Let us start by introducing the Black and Scholes model driven by the Rosenblatt 
process. As usually, we will consider two assets: a safe investment satisfying 



TsBsds, 



(15) 



where r is a differentiable deterministic function and a risky asset with price dynamic following 
the stochastic equation 



St = So+ / ttsSgds + <y I SgdZs 
Jo Jo 



(16) 



where Z is a Rosenblatt process, cr > and a is a differentiable deterministic function. The 
integral with respect to Z is understood here in a pathwise sense. Since the trajectories of 
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the Rosenblatt process are enough regular (in fact, they are Holder continuous of order 6 < H 
and H > ^) it is possible to consider pathwise integrals with respect to it and to solve some 
stochastic equations in the pathwise sense. We refer, among others, to [13], [27] or [T2]. In 
particular, the solution of (fT6]) is given by (see [H], [27], [13] 

St = Soe-^o ^sds+aZt ^ t g [0, T] . (17) 

Clearly the solution of (fT5|) is 

Bt = Boe-^o^^'i'^ tG[0,T]. (18) 

Moreover, we will assume in the sequel that the interest rates a and r are deterministic bounded 
functions. 



Let us describe now the binary market model with Rosenblatt influence. The two assets 
are traded now at successive times periods ti < t2 < • • • < ••• and their dynamics are given by 

Bn = {l + rn)Bn-l 

and 

Sn = {an + (1 + Xn))Sn-i; 

That means that Bn and Sn represent the prices of the bond and of the stock in the period 
between tn and tn+i and r„, and a„ are the interest rates valuable in this period. The stochastic 
process X is binary, that is, given Xn-i it can takes at time n two possible values denoted by 
d^C'down") and Un ("up"). The binary market excludes arbitrage opportunities if for every n 
it holds that (see |22j) 

dn < Vn - an < Un- (19) 

In the following we will choose our binary model to be determined by 

Xn = AZ^ (20) 

where is defined by (jlip and 



— riL and a„ = — an (21) 



where a and r are the interest rates appearing in (jlSp and ([16 



We have 

Proposition 3 The binary market model with X a and r given by I12(J\) and I121\) converges as 
N ^ oo to the Black and Scholes model given by U5\) and il6\) . 



10 



Proof: Let us consider the jump 



and the quadratic variation 



[z% = Y^i^z^ 



N\2 
s ) ■ 



s<t 



We win show that the process [Z^] converges in L-^([0,T] x Q) to zero. Then the conclusion 
will follow exactly as in |22], proof of Lemma 1. 

We have, since the jumps are at times jj^, k integer, 



and then 



This implies that 



E lAZf 1^ < E 



E[Z^]t < Nt 



- Z,i 



2 1 

< 



tN 



1~2H 



rp 

[ [Z%ds<c{T)N 
Jo 



1-2H 



which goes to as goes to oo. 



The next step is to show that the market admits arbitrage opportunities. Clearly, we 



have 



and we will take the random variables ^ to be binary, that is 

P(e = 1) = = -1) = ^ 

for every i > 1. We can write, by isolating the part involving 

Xn = /n-l(6> • • • >Cn-l) + '^n5n-l('^l, • • • , (,n-l] 

where for every n>2 

n — 1 



J. ^ J. / 

N f N f / n \ ^ 



M, t> I I dvdu XiXj 
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and 

n-l 



Then obviously 



u,v] dvdu 1 Xi 



Un = fn-l{^l, - ■ ■ +9n-l{Cl,---,^n-l) (22) 

and 

dn = fn-l{£,l, ■ ■ ■ ,?n-l) " 5'n-l(Cl> • • • :Cn-l)- (23) 

The last result of our paper is the following. 
Proposition 4 The binary market model with Ii21\) and Ii20\) admits arbitrage. 



Proof: Throughout this proof, we will denote by c{H) a generic constant depending only on 
H. Let us show now that the condition (jl9p fails for some n > 2. We will actually prove that 
the sequence 

1, . . . , 1) - gn-l{l, 1, • • • 1) ^n-*oo OO (24) 



and then clearly (|T9|1 does not hold because r„ and an are assumed to be bounded. 
We have 

!,...,!) = aN J2 ( r. (^(^'"'^)-^(^'^'^))^"^^ 

= E (/!/*. ('^(^■''•")-''(^.«'''))") 

71—1 / „J_ / / 

N f N f ^ / n \ ^ f n — 1 11, , 
V ] — F \ — T^,u, V I I dvdu 



i=l \ N N 

aN{A-B) 
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Using the expression of the kernel F and ([U]), the term A can be minorized as follows 

n—1 n— 1 n 



A = c{H) / / (a-u)^ -2{a-v)^ -2a^^ ~^U2-^ V2-^ dadvdu 

Jo Jo 



n-1 
N 



. , 2//' n — 1 n — 1 n 



We majorize now the term B. We can write 



n — 1 „_L 



^ / / . .u'_3, ._H-'_3 2H'-1 i-H' 



2 n2 ^ t;2 ^ dadvdu 



^ — , f N r N r N rji 3 I 

« £ <^mE/_. /_,/._. («-")" 

/ 72 — 1 _i_ _L ZL 



i=l AT iV JV 

1 = 1 — — 



< <H) - 



NJ ^\\Nl \ N 1 I iV2H'-2 

i=l 

n\2H'-l 1 /^_p2H'-2 



From the above computations we obtain that . . . , 1) converges to oo as n ^ oo because 

2H' — 1 > AH' — 3. Let us estimate the term (7„_i(l, . . . , 1). 
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£/n-l(l, •••,!) 

i=i \-^— -^^r / 

/ 'M— 1 _L — — 

, N / n \ -1 ^ ^ f N f N f N \ Tjl 1 TTl ,n — 1 ^H■'3, NW'f 

< 2cjiVc(i/) (^— j ^ , , {— u)^ -2{a-v)" -2dadvdu 

i=i -^hr "^TT" 

< 2aNciH) (-) (^^j 5: y^_^ u.-- - -H- - '.dvdu 

^ ^ %=!•'—•' -N~ 

< 2<TjV<;(ff)^^^^o(„«'-|) 

because the integral Jq ^ u2~^ {^^ — u)^ ~2du is equal to — i/', H' — \). We then obtain 
(IMI). ■ 



Comments: i) A concrete arbitrage opportunity can be easily described. For exam- 
ple, suppose that a > r. Ar a certain time no we have dng > because of (fM|) : suppose that 
the stock price was increasing up to this time no- Then, buy M stocks and your wealth at 
time no + 1 will be positive since MSuq+i > MSuq- 

ii) We cannot expect to have no-arbitrage when G ^ as in the fractional Brownian 
motion case because now the limit process at = ^ is not necessarily a martingale. 
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5 Appendix: Monte Carlo Simulation 



This simulation study is intended to show our proposed model using simulated data. We 
simulate data using different values for H. Below, we describe the procedure used in generating 
the data to be used in the simulation study. 

We have implemented this simulation on a standard personal computing platform (PC), 
and have observed that it performs very well using simulated data as can be seen from the 
simulated data in the figures 1 and 2 below. Despite the apparent algebraic complexity of 
the equations (jllh . the problem poses no difficulty for standard symbolic algebra packages. 
Using Matlab's simulations and algebra capabilities yielded the best computing times. In our 
implementation, which performs an iteration of the algorithm from i = to i = n. Figure [1] 
shows the histogram for a fixed time t for the marginal density of the Rosenblatt process. We 
can see the skew structure of the distribution. 






-0.5 0.5 



2.5 3 3.S 



Figure 1: Marginal distributions for H = 0.5; H = 0.8 and H = 0.9. 

Figure [2] shows the some paths of the discretization for the Rosenblatt process. We use 
the values for the parameter H (H = 0.8 and H = 0.9). 





Figure 2: Simulations of the Trajectory for the Rosenblatt process with H = 0.8 and H = 0.9. 
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